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Abstract 

We investigate the generalization of the Costas property in 3 or more dimensions, and we seek an appropriate 
definition; the 2 main complications are a) that the number of "dots" this multidimensional structure should 
have is not obvious, and b) that the notion of the multidimensional permutation needs some clarification. After 
proposing various alternatives for the generalization of the definition of the Costas property, based on the 
definitions of the Costas property in 1 or 2 dimensions, we also offer some construction methods, the main one 
of which is based on the idea of reshaping Costas arrays into higher-dimensional entities. 

1 Introduction 

Costas arrays originated in engineering as a frequency hopping pattern that optimizes the performance of radars 
and sonars l3]|4]; being a singular combinatorial object, however, they have been lately the focus of intensive 
study by mathematicians as well, and have thus started leading a second independent "life" in the mathematical 
literature [711511^ [1111114) . In this work, we increase the level of mathematical abstraction by investigating analogs 
of Costas arrays in 3 or more dimensions: several challenges lie ahead, as we first need to give a satisfactory 
definition of the Costas property in higher dimensions, and also provide some algorithms to construct such 
Costas cubes or "hypercubes" in general. 

We will begin by stating the definition of Costas arrays [5j in such a way as to exhibit their direct relation 
to Golomb rulers [16], which will consequently be construed as the 1-dimensional analog of Costas arrays. 
Subsequently, we will discuss several slightly different candidate definitions in higher dimensions, and we will 
attempt to generalize the generation methods available to us. We shall see, in particular, that it is possible to 
generate sparse Costas hypercubes systematically, using various construction methods, but if we want to know 
what the densest Costas hypercube can be, things become less clear. 

2 The definition of Costas property in 1 and 2 dimensions 

We use the usual notation that [n] := {1,2, . . . ,n}, n G N, and that [n] — 1 := {0, 1, ... ,n — 1}, n G N*. 

Definition 1 (Costas property in 1 dimension). Let g : [N] — » [N'], N,N' £ N, be a finite strictly increasing 
sequence, such that g{l) = 1, g{N) = N'; g will satisfy the Costas property in one dimension iff g{i) — g{i + k) — 
g{j) — g{j + k) i = j, i,j, i + k,j + k G [A'^]. Equivalently, if / : [A*''] — > {0, 1} is a sequence such that f{i) = 1 -^^^ 
3j G [N] : i — g{j), and f{i) = whenever i ^ [N'], then the autocorrelation Af{k) = X]ig[jv'] + k & Z 

satisfies Af{k) < 1, k ^ 0. The sequence g is known a Golomb ruler [16j . 
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Definition 2 (Costas property in 2 dimensions). Let ] -.l? ^ {0, 1} be a binary sequence in 2 dimensions, such 
tliat it is equal to 1 at a finite number of points only; without loss of generality we may assume that f{i,j) ~ 
when i ^ [M] or j ^ [A'"] for some M,N £ N. Defining the autocorrelation as Af{k,l) — j^j^ f{i, j)f(i + k, j -\- 
l), k,l £ Z, f will have the Costas property in 2 dimensions iff Af{k, I) < 1, k,l ^ 0. 

Definition 3 (Costas squares). Let / have the Costas property in 2 dimensions and let M,N, appearing in 
Definition [2] be the least possible; if M = N = n, f will be called a Costas square of side length n. 

Definition 4 (Costas arrays). Let / be a Costas square, and let it, in addition, have the structure of a 
permutation matrix with only one element equal to 1 per row and column: in other words, let there be a 
bijective sequence g : [n] —> [n] (namely a permutation of order n) such that f{i,j) = 1 <^ i = g{j),i,j G [n]. In 
that case, / will be called a Costas array, and g a Costas permutation. 

Remark 1. 

• It follows immediately that, if g is a Costas permutation of order n, no two vectors in the collection 
{(* ~ Ji S'(*) ~ 9ij)) • 1 ^ * < j ^ '^j can be equal, and that no vector in this collection has a coordinate 
equal to 0. 

• It is customary to denote Os in the array / by blanks and Is by dots. 

• Note that the Costas property, as stated in Definition [J] is (much) more general than usual, as it can be 
satisfied by (many) more 2-dimensional sequences than Costas arrays themselves, which appear only as 
quite restricted special cases (see Definition |4ll . But this generality will actually prove to be helpful when 
we seek a higher-dimensional analog of Costas arrays. 

3 Costas hypercubes 

It is straightforward to generalize the Costas property in higher dimensions, to the extent that we risk (re)stating 
the obvious below. What is harder is to find the exact higher-dimensional analog of the special case, the "Costas 
array" . 

3.1 The Costas property in higher dimensions 

Definition 5 (Costas hyper-rectangles). Let m G N, consider a sequence / : {0, 1}, and suppose further 

that f{i) =0, i ^ [N],N G N"*, where i = {ii,...,im), N = {Ni,...,Nm), [N] = [iVi] x ... x [N^], and 
the vector A'^ has the smallest possible entries (for the given sequence /). Let the autocorrelation of / be 
^/(fe) = Eiez". /(«)./"(« + ke Z™. Then, / will be a Costas hyper-rectangle iff Vfc G Z"" - {0}, Af{k) < 1. 
Definition 6 (Costas hypercubes). If / is a Costas hyper-rectangle so that A^'i = . . . = Nm = n G N, it is called 
a Costas hypercube; if m = 2, it is a Costas square (see Definition |3]), whereas if m = 3 it will be called a Costas 
cube. 

Now it is time to seek the higher-dimensional analog of a permutation. 

3.2 Vector permutations 

It is quite simple to generalize permutations in higher dimensions, if the number of dimensions is even: 

Definition 7 (Permutation Costas hypercube). Let m — 2s, s G N, and let g : [n]" [n]" be a bijection, 
that is a permutation on vectors in general. Let / : — > {0, 1} be a sequence such that f{i) = 1 iff 
{is+i, . ■ . , i2s) = g{ii, . . . ,is), (ii, . . . , is) G [nY , and such that it has the Costas property, as defined in Definition 
[S] then, / will be called a permutation Costas hypercube in m dimensions with side length n. 

Remark 2. 

• Permutation Costas hyper-rectangles are defined by the obvious extension of the above definition. 

• The fact that we chose the first s dimensions to form the domain of g and the last s its range does not 
affect generality: if / is a Costas hypercube, then any /' resulting by a random permutation of the order of 
the dimensions is also a Costas hypercube; this constitutes a generalization of the invariance of the Costas 
property in 2 dimensions under transposition. 

• No 2 vectors in the collection {{i — j, g{i) ~ g{j)) '. ,i 7^ j} can be equal; however, they may have 
coordinates equal to 0. 
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Table 1: A hypercube with m = 3, n — A constructed by permuting all pairs of integers 

Example 1. As a specific example, let us use m = 4, n = 3. Then, the hypercube with f{i) = 1 iff i is one of 
the row vectors of Table [32] is a permutation Costas hypercube; this can be checked by a) verifying that the 2 
first columns contain all vectors with integer coordinates between 1 and 3, as do the 2 last columns, and b) by 

finding all possible 11 =36 distance vectors and observing they are indeed distinct. Observe also that this 

hypercube has in total 9 = 3^ nonzero elements out of 81 — 3*; in general, permutation Costas hypercubes have 
nonzero elements out of n^" . 

Remark 3. When m = 2s + 1, s £ N, it is clearly impossible to define a permutation as we did above; the 
best we can aim for is to find an injective function g : [n]" — > n"'^^ , so that f{i) = 1 iff {is+i, ■ ■ ■ ,i2s+i) = 
g{ii, ■ ■ ■ , is), {ii, . . . ,is) £ [n]" ■ Obviously, any Costas hypercube in m — 1 = 2s dimensions can be extended to 
m dimensions, by adding one more coordinate to the position vectors and assigning values randomly to it. This 
can be done in n"^^ ways, as there are n" vectors and in each one the last (new) coordinate can now take any 
value in [n]: the subvectors formed by the first 2s coordinates of the distance vectors remain distinct, so the last 
coordinate is effectively unused (and therefore purely ornamental). It is of course possible to turn hypercubes 
in 2s dimensions that do not have the Costas property into Costas hypercubes by a judicious choice of the new 
coordinate's values. 

3.3 Strict Costas hypercubes 

Vectors between dots in a permutation Costas hypercube are allowed to have coordinates, in contrast to 
ordinary 2-dimensional Costas arrays where this does not happen (see Remark We could potentially restrict 
the definition of a Costas hypercube to exclude such a possibility, although such a restriction is extremely severe: 
invoking the Pigeonhole Principle, if a hypercube of side length n has n + 1 dots, then, fixing any coordinate, 
at least 2 of those have position vectors with the same value for the chosen coordinate, hence the corresponding 
distance vector has a there. Therefore, our requirement limits the number of coordinates/dimensions to at most 
n, implying that one of the coordinates (of an element of the hypercube equal to 1) determines unambiguously 
the rest of them. 

Definition 8 (Strict Costas hypercubes). Let m, n G N and let / : {0, 1} be such that f{i) = 1 iff 

i G {ii, . . . ,in}, where ij £ M'", j £ [m], and where all vectors in the family {ij — ifejl < j < fc < n} are 
distinct and have no coordinates equal to 0; then, / will be called a strict Costas hypercube. 

We proceed to give 2 explicit construction methods for strict Costas hypercubes. 

Theorem 1 (Lifted Costas hypercubes). Let m, n G N and let f : ^ {0)1}- Choose m — 1 permu- 
tations Qi, i £ [m — 1] (not necessarily distinct), so that g\ is in addition Costas, and set f[i) = 1 iff 
i G {(ii gi{j)y • • • 1 9m-i{j)) '■ j = 1, ■ ■ ■ , n}. Then, f is a strict Costas hypercube. 

Proof. The result is practically obvious: a typical distance vector is {j — k, gi{j) — gi (fc), . . . , pm-i (j) — ffm-i(fe)), 
where 1 < j < k < n, and as all gs are permutations, no coordinate is equal to 0. Further, the first 2 coordinates 
of the vectors above are all distinct, as gi is a Costas permutation. □ 

Theorem 2 (Toeplitz-Costas hypercubes). Let m, n G N, n < m, and consider the column vector u = 
(1, . . . , n)', the vertical shift operator S{xi, . . . , Xn)' ~ {xn, xi, . . . ,a;„_i)', and the array A = [u Su . . . S"^~^u] . 
Let f : Z'" {0, 1} and set f(i) = Q iff i is one of the rows of A. Then, f is a strict Costas hypercube and A 
its corresponding Toeplitz array. 
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1 4 3 2 1 

2 14 3 2 

3 2 14 3 

4 3 2 1 4 



Table 2: A Toeplitz construction of a Costas hypercube with m = 5, n = 4 



Proof. Every column of A contains each integer in [n] exactly once, hence no distance vector has a coordinate 
equal to 0. Further, given the difference between 2 rows, the 2 rows can be uniquely determined. Indeed, if r; 
and Tj denote the rows i and j of A, respectively, with i < j, then — rj is a vector with the first i coordinates 
equal to i — j, the following j — i coordinates equal to n + i — j, and the remaining coordinates equal to i — j; 
hence, i and j is uniquely determined by inspection. □ 



Example 2. Table 1331 shows the corresponding array A of Costas- Toeplitz hypercube with n — A and m = 5. 

The former method has the drawback that it requires that Costas permutations of order n be known, while 
the latter has the drawback that m > n, so in general the number of dimensions m needs to be very high; on 
the other hand, it requires no Costas permutations of order n to be known, so it can be used for values of n 
such as 32 and 33, where no Costas arrays are known yet. In any case, strict Costas hypercubes have lots of 
"empty space", namely an extremely low density of dots (they only have n dots), and therefore they tend not 
to be very interesting. 

4 The main construction method 

It is natural to ask whether Costas arrays (in 2 dimensions) can somehow be manipulated (essentially reshaped) 
to produce Costas hypercubes, and hopefully permutation Costas hypercubes. We prove below that this is 
possible, and provide several variants of this construction. 

4.1 Reshaping 

We formulate and prove below a general result about constructing a Costas hyper-rectangle out of a Costas 
square. In the special case where the Costas square is a Costas array, and the hyper-rectangle a hypercube, it 
turns out the hypercube is a permutation Costas hypercube. 

m 

Theorem 3 (Reshaping). Let m,n £ W , n = Y\ 'n-i, n-t > 1, i G [m], and let g he a Costas permutation of order 

m m 

n, hut following the convention that g : [n] — 1 ^ [n] — 1. Expand i — ^ ^ Vj (i) ni, so that i gets mapped 

j = l 1=3 + 1 

hijectively to V{i) — {vi{i), . . . ,Vm{i)), where Vj £ [nj] — 1, j € [m]; similarly, g{i) gets mapped hijectively to 
V{g{i)) = (vi{g{i)), . . . , (<?(«)))• Then, the hyper-rectangle of side length m in dimension i and i + m, i £ [m], 
whose dots (n m total) lie at the points {V{i), V(g{i))) :— {vi{i), . . . , Vm{i), vi{g{i)), . . . , Vm{g{i))), i G ["■] ^ 1; 
is actually a permutation Costas hyper-rectangle. 

Proof. 

1. Choose 2 values for i, say ii and i^; the corresponding distance vector is: 

{V{ii) ~V{i2),V{g{t^)) -V{g{i2))) = 

= - -^1(42), . . . ,Um(ji) - v,n{i2),vi{g{ii)) - vi{g{i2)), . . .,Vm{g{ii)) - v^{g{i2))) 

We need to show that all of these vectors are distinct. In other words, we need to show: 

(V(ii) - V{i2),V{g{i^)) - V{g{i2))) = {V{iz) - V {i^) ,V {g{iz)) - ^(5(^4))) ^ ii = ^2, h = i4 

m m 

2. Extend by {vi, . . . ,Vm) = '^i Y\ class of vectors where \vj\ < uj, j G [m]. It 

i=i i=j+i 

follows that V~^{V{ii) — V{i2)) — ii — 12, as V{ii) — V{i2) falls within this class of vectors. 
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3. Putting things together: 



iVin) - V{i2),V{g{ii)) - V{g{i2))) = {V{h) - V{U),V{g{is)) - V{g{U))) ^ 

V-'{V{h) - V{i2),V{g{h)) - V{g{i2))) = V-\V{ia) - V{u),V{g{h)) ~ V{g{u))) ^ 
{V~\V{h) ~ V{i2)),V''{V{g{H)) - V{g{i2)))) = {V-\V{h) - V{u)),V''{V{g{i,)) - V{g{u)))) ^ 
{ii-i2,g(ii)~g(i2)) = {i3-ii,g{i-A)-g{i4,)) <4> 11-12 = h-ii, g{ii)-g{i2) = g{i3)-g{u) => ii = i-i, 12 = ii 

In the last step we used the fact that g is Costas. 

Further, observe that the left half row vectors so constructed are the expansions of i G [n] — 1, while the 
right half vectors are the expansions of i £ [n] — 1; the fact that every i £ [n] — 1 gets expanded 

exactly once and that g is a permutation guarantees that the hyper-rectangle produced is a permutation 
one. This completes the proof. 

□ 

The hypercube of even dimension is just a special case: 
Corollary 1 (Costas hypercube of even dimension). Let m,n £ N*, and let <; be a Costas permutation of order 

m 

n™, but following the convention that g : [n]™ — 1 ^ [n]™ — 1. Expand i = Vj{i)ri^~\ so that i gets mapped 

i=i 

bijectively to V{i) = (i'i(i), • • • ,Vm{i)), where Vj G [n] — 1, j G [m]; similarly, g(i) gets mapped bijectively to 
V{g{i)) — {vi{g{i)), . . . ,Vm(g{i)))- Then, the hypercube of side length n whose dots (n™ in total) lie at the 
points (V(i), V{g{i))) :— {vi{i), . . . , Vm{i),vi{g{i)), . . . , Vm{g{i))), i € [n™] — 1, is actually a permutation Costas 
hypercube. 

Remark 4. Notice that these hypercubes have n"* dots, which is the square root of the n^"* total positions 
available in the hypercube, just like in Costas arrays, where there are n dots among the available positions. 

Ifere is an attempt to construct approximate Costas hypercubes of odd dimension in the special case where 
the side length n is a perfect square, using Theorem [S] by first constructing a hyper-rectangle as an intermediate 
step: 

Heuristic 1 (Costas hypercube of odd dimension). Let m,n G N*, where ^/n G N, and let g be a Costas 
permutation of order n"^^/n, but following the convention that g : [n]"* x [-y/rll — 1 ^ x ^ 1- Expand 

m 

i = vo{i)n"^ + y ^ Vj (i)n'"~"' , so that i gets mapped bijectively to V{i) = (wo(i), Vi{i), . . . , Vm{i)), where Vj G [n] — 
i=i 

1, j G [m] andvo G [y'n] — 1; similarly, g{i) gets mapped bijectively to V{g{i)) — (yo{g{i)),vi{g[i)), . . . ,Vm{g{i)))- 
This process forms a Costas hyper-rectangle in 2m + 2 dimensions, whose side length in 2m dimensions is 
n and in the remaining 2 dimensions ^/n. Now, replace the coordinate pair {vo{i),vo(g{i))) in the coor- 
dinate vector of each dot by the single coordinate ^/nvo{g{i)) -\- vo{i), i G [n"^^/n\ — 1, which takes values 
in the range [n] — 1. Then, the hypercube of side length n whose dots (ri^^/n in total) lie at the points 
(x/nfo ((?(*)) vo{i),vi{i), . . . ,Vm{i),vi{g{i)), . . . ,Vm{g{i))), i G [n^y^l ^ 1; usually a good approximation 
of a Costas hypercube, and the removal of a few dots turns it into a Costas hypercube. 

Remark 5. The reason why this heuristic often fails to produce a Costas hypercube is that different pairs 
of coordinates (in difference vectors) can collapse to the same value: for example, assume n = 25 -^^ ^/n = 5 
and consider the pairs (—3,0) and (2,-1); they get mapped to —3 + 5 • = —3 and 2 — 5 = —3. In the 
context of the proof of Theorem |3l this is equivalent to saying that V{ii) — V{i2) is not always the same as 
sign(ii — i2)V{\ii — 12]). However, simulations show that the damage this does to the Costas property is usually 
small: tests with Costas arrays of side length n"^ < 200 showed that systematically over 95% of the difference 
vectors among the dots are distinct. 

The construction methods above can be significantly extended if we use Costas squares instead of Costas 
arrays as the starting point. The key observation (through the proof of Theorem |3} is that the permutation 
property of the original Costas array is not responsible for the Costas property of the hyper-rectangle produced, 
but rather for its permutation property alone (and in the case of Heuristic [T] it does not even achieve that). 
Therefore, if we are not interested in obtaining a permutation Costas hyper-rectangle as the final product, or if 
a suitable sized Costas array is not available, we may as well start with a Costas square. 

A special type of Costas squares that proves very helpful in practice is smaller Costas arrays. Consider a 
Costas array of order n' G N* and let n > n'; then, this Costas array can be turned, by the addition oi n — n' 
blank rows and columns at the sides of the array, into a Costas square of size n. Note that such Costas squares 
are generated by incomplete permutations. We generalize this notion in the following definition: 



5 



Definition 9 (Incomplete Costas array). A Costas square with the property that there is at most one dot 
per row and column will be called an incomplete Costas array. A (Costas) hyper-rectangle/hypercube of even 
dimension with the property that there are no 2 dots whose position vectors have the same left half or right half 
part will be called an incomplete (Costas) hyper-rectangle/hypercube. 

Corollary 2 (Constructions out of Costas squares). 

• A construction according to Theorem[3]starting with a (incomplete) Costas square results to a (incomplete) 
Costas hyper-rectangle. 

• A construction according to Corollary[T]starting with a (incomplete) Costas square results to a (incomplete) 
Costas hypercube. 

• A construction according to Heuristic [T] starting with a (incomplete) Costas square results to a hypercube 
that very nearly has the Costas property and usually can be turned into a Costas hypercube through the 
removal of a few dots. 

4.2 Construction examples 

We give 3 examples: the first is the construction of a Costas hypercube with n = 5, m = 4 out of a Costas 
array of order 25 (using Corollary [1} ; the second is the construction of a Costas hypercube with n = 9, m = 3 
out of a Costas array of order 27 (using Heuristic [1} ; and the third is the construction of an incomplete Costas 
hypercube with m = 5, n = 4 (using Heuristic [T]) , starting with a Costas array of order 31 and extending it into 
an incomplete Costas array of order 32. 

Example 3. Consider the permutation of order 25 appearing on Table |3] (left). Applying Corollary [TJ we get 
a Costas hypercube with m = 4, n = 5, whose dot positions appear on Table [3] (right). Observe that this is 
a permutation Costas hypercube, as Corollary [T] states: every vector (i,j), i,j G [5] — 1 appears in the left 2 
columns in exactly 1 row and in the right 2 columns also in exactly one row. 

Example 4. Consider the Costas permutation of order 27 appearing on Table |4] (left). Applying Heuristic [T] 
we first get a Costas hyper-rectangle in 4 dimensions of side lengths 9 and 3, whose dots lie at the points shown 
in Table m (center) . Subsequently, we combine the middle columns that have values in the range {0, 1,2} into a 
single column with values in the range [9] — 1, as described in Heuristic[T] for example, the middle 2 coordinates 
(1,2) in row 12 become 1-1-2-3 = 7. The result appears in Table (right) . A check of the Costas property 
shows that this time we get lucky and that the hypercube of m = 3, n = 9 we have created is Costas. It is 
obviously not a permutation Costas hypercube, as this term is meaningless in odd dimensions. 

Example 5. Consider the Costas permutation of order 31 appearing on Table [S] (left), and consider it as an 
incomplete Costas permutation of order 32 = 2^ = 4^\/4- Applying Heuristic [TJ we first get a Costas hyper- 
rectangle in 6 dimensions of side lengths 4 and 2, whose dots lie at the points shown in Table [5] (center) . 
Subsequently, we combine the middle columns that have values in the range {0, 1} into a single column with 
values in the range [4] — 1, as described in Heuristic [1] The result appears in Table [5] (right). A check of the 
Costas property shows that this time we get lucky and that the hypercube of m = 5, n = 4 we have created is 
Costas. It is obviously not a permutation Costas hypercube, as this term is meaningless in odd dimensions. It 
also has just less than \/4^ = 32 dots, one less than that to be exact. 

4.3 Applicability of reshaping 

The construction method we described reshapes a Costas array into a high-dimensional Costas hypercube. As 
the order n™ is bound to be quite big, the only Costas arrays that will normally be available are Golomb and 
Welch constructions [5] . For practical purposes, we can seek Costas arrays of large orders in databases, such as 
the database created by Dr. J. K. Beard containing all known Costas arrays up to the order 200 ,2 . But more 
generally, in order to figure out whether a 2m-dimensional cube of side length n exist, we will need to check 
whether one of the following equations holds: 

• -f 1 = p (Wi construction possible) 

• n™ -I- 2 = (G2 construction possible; if fc = 1, 11^2 construction possible, too) 

• -f 3 = (G3 construction possible) 

where in all cases p is a prime. Other variants of the Golomb and the Welch construction do not occur 
systematically, so we do not investigate them. We now look briefiy into each one of the above equations. As a 
general comment, the solution of these equations falls under the scope of Diophantine Analysis, and it appears 
that many conjectures can be formulated, but few facts have actually been proved. 
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Table 3: The conversion of a Costas array of order 25 into a Costas hypercube with m = 4, n = 5: the permutation 
(left), and the final Costas hypercube (right) 
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Table 4: The eonversion of a Costas array of order 27 into a Costas hypercube with m = 3, n = 9: the permutation 
(left), the intermediate hyper-rectangle (center), and the final Costas hypercube (right) 
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Table 5: The conversion of a Costas array of order 31 into a Costas hypercube with m = 5, n = 4, treating the 
Costas array as an incomplete Costas array of order 32: the (incomplete) permutation (left), the intermediate 
hyper-rectanglc (center), and the final Costas hypercube (right) 
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4.3.1 + 1 = p 

Assume that m — m\ra2, where mi is odd; then nJ^^ + l|n™ + 1, hence it cannot be a prime. It foUows that 
m = 2 , for some fc, whence p — n +1. Further, n must necessarily be even, whence n = 21 and p = 2 I + 1. 
For Z = 1, we obtain p = 2^ +1, the celebrated Fermat primes, for which a lot is known; in particular, it is 
conjectured that the only such primes correspond to 1 < fc < 4, and lead to hypercubes of side length 2 in 2^* 
dimensions. 

Let us now fix fe = 1, namely seek primes of the form p = 1 + n?. It is conjectured there are infinitely many 
such primes [11], and it is easy to find several examples: 17 = 4^ + 1, 37 = 6^ + 1, 101 = 10^ + 1, 197 = 
14^ + 1, 257 =16^ + 1 etc. 

Solutions for higher k can also be found. For example, for 1 < n < 20 and A; = 2, n* + 1 is prime 
for n = 2, 4, 6, 16, 20. Those solutions can actually build several hypercubes, namely either of side in 4 
dimensions, or of side n in 8 dimensions, in all cases with dots. 

4.3.2 n™ + 2=p'= 

Here necessarily n must be odd, and again solutions can be found: 3^ + 2 = 11, 9^ + 2 = 83, 15^ + 2 = 227; 
3^ + 2 = 29, 5^ + 2 = 127; 5^ + 2 = 3^ etc. It is not known whether infinitely many solutions exist. This 
equation appears to be a generalization of Catalan's Equation [13j . 



4.3.3 n" + 3 = p*^ 

Here n must be even, and solutions also exist: 2* + 3 = 
It is not known whether infinitely many solutions exist. 
Catalan's Equation jl3i . 



19, 2^ + 3 = 67, 10^ + 3 = 103, 14^ + 3 = 199 etc. 
This equation also appears to be a generalization of 



4.4 Older generalization attempts 

Extensions of the Golomb and the Welch constructions to 3 dimensions were investigated in the past [6] , but the 
objective then was slightly different: the cubes were so constructed that all the 2-dimensional "slices" along some 
of their directions be Costas arrays, or at least almost Costas arrays in a certain sense. In other words, although 
the construction was 3-dimensional, the Costas property was still investigated in 2 dimensions. It turns out that 
many "reasonable" generalizations yield cubes with with 2-dimensional (almost) Costas slices. For example, the 
dots in the cube are placed at the points whose coordinates are the solutions k), i,j, k £ [q — 1], of one of 
the following equations: 

1. a*"""^ + b''^'" + c^^'^ = 0, where a, 6, c are primitive roots of the field F((j'), for some q > 2 power of a prime, 
and X, y, z fixed in the range 0, . . . , g — 2; 

2. a*"""^ + = k mod p, where p prime, a, 6 primitive roots of the field F(p), and x^y fixed in the range 
0,...,g-2; 

3. a*"'"-' = k mod p, where p prime, and a primitive root of the field F(p); 

4. a*"'"'* = jk mod p, where p prime, a primitive root of the field F(p), and d fixed in the range 0, . . . , g — 2. 

Note that the first one can be expanded in an obvious way into arbitrarily many dimensions. Unfortunately, 
computer simulations show that none of these constructions yields a Costas cube, so they are not suitable for 
our purposes. 

5 An extension of the Welch construction 

5.1 The original method and its non-extendability for Costas arrays 

The Welch construction method for Costas arrays [5] |7] stipulates that, if p is a prime, g a primitive root [T] 
of F(p) and c G [p — 1] — 1 a fixed parameter, then the function = g'~^^'^ mod p, i G [p — 1] is actually a 
Costas permutation on [p — 1]. Contrary to the Golomb construction [5] [7], though, which works in all finite 
fields, namely with p'" elements where p is a prime and m £ N*, the Welch method is not applicable when 
m > 1. 

The reason for that is simple: when m > 1, the elements of the field are no longer represented by integers, 
but rather by polynomials of degree m — 1 in a (for all practical purposes) "dummy" variable, say x, and with 
coefficients in [p] — 1, while addition and multiplication are no longer defined modulo an integer, but rather 
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modulo a monic irreducible polynomial P{x) of degree m [T] [5]. It follows that the function / of the Welch 
construction is now / : [g — 1] — > F*(g) where f{i) = g^~^'^'^ mod P{x), i £ [q — 1], with c £ [g — 1] — 1 and g 
a primitive root of ¥{q) (it can be shown that the multiplicative subgroup F*(g) of the field F(g) is still cyclic 
[T][5]); therefore, f{i) is a polynomial, while i is an integer, and the whole construction is (at first sight, at least) 
meaningless, as we need / to produce integer values! 

Perhaps a more correct way to think of the variable x is that it represents an algebraic element of F(p) of 
order m; that is, it is a root of the polynomial P{x) with coefficients in F(p) and of degree m, while it is the 
root of no other (non-zero) polynomial with coefficients in F(p) and of degree less than m [T]. 

Although this construction fails to yield a Costas array, it still produces Costas hyper-rectangles, as we are 
about to see. 

5.2 Construction of hyper-rectangles 

The field F(p'") can be construed to be a vector space over the field F(p) in 2 ways, depending on whether we 
consider its elements to be polynomials or m-tuples: 

Definition 10. The field F(p'"), where p prime and m G N*, when viewed as a vector space over the field F(p), 
will be denoted by F(p) . Let Vp{p,m) denote the vector space of polynomials of degree m—1 over the field F(p) . 
Then, Vp{p, m) and F(p)™ are isomorphic vector spaces, as they have the same finite dimension m and they are 
over the same field p^, under the isomorphism denoted by whereby T{am-ix"^~^ + am-2X^~^ -I- ... -I- ao) = 
(am,-i,am_2, . . . ,ao), G F(p), i G [m] — 1. 

Theorem 4 (Welch hyper-rectangles and hypercubes). Let p he a prime, m £ W , g a primitive root of¥{q) 
where q = p™", and c G [g — 1] — 1. We choose Vp{p,m) as our representation of¥{q). Then: 

• The function / : [g — 1] ^ F*(g), where f{i) — jr((;'^^+'^ mod -P(a;)) , i G [g — 1], P{x) an irreducible 
polynomial over F(p) of degree m, is a permutation over F*(g). 

• The hyper-rectangle in m + 1 dimensions with side length g — 1 in the first dimension and p in the others, 
whose dots lie at the positions with coordinates {{i, f{i))\i G [g — 1]}, has the Costas property. 

• The hypercube in 2m dimensions with side length p, whose dots lie at the positions with coordinates 
{{V{i), f{i))\i G [g — 1]}, has the Costas property; V is the familiar mapping from Corollary \1\ with 
n = p in the present case. 

Proof. The proof really consists of putting together bits we have already proved. 

• ,f is a permutation over F*(g) because g is taken to be a primitive root of F(g). 

• That the family of vectors {{i, f{i))\i G [g— 1]} has the Costas property follows from a verbatim repetition 
of the classical argument for m — 1 |5! "7]. Let us see it here in detail: consider the 4 integers ii, 12, 
is = ii + k, and 14 = 12 + fc such that ij G [g — 1], j G [4], fc G N, ii < 22, and let (ii — 12, f{ii) — f{i2)) ~ 
(is — 14, /(is) — f{ii)) <^ ^ /(ia) = /(is) — fiu). This can be written as: 



Hence, ii — is, i2 = ii and the proof is complete. 

• We need to show that, given that the family {{i, f{i))\i G [g — 1]} has the Costas property, the family 
{(V(i), f{i))\i G [g— 1]} has it too; but this is a verbatim repetition of the argument presented in the proof 
of Theorem |3] 



• The hypercubes constructed above have g — 1 dots out of g^ possible dot positions, and thus follow 
approximately the square root rule we saw earlier for the density. 




(g*" - l)(p'i -p'2) = mod P(a;) <^ 3''' - 1 = mod P(x) V g'^ - g'" = mod P(j:) ^ 
g'' = 1 mod P{x), as < ii < 12 < g, which makes the second condition always false 



fc = 0, as < fc < g - 2. 



□ 



Remark 6. 
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• They are "almost" permutation hypercubes, except that the zero vectors are missing. If we set /(O) = 
and add a dot at the position (1/(0), /(O)) — (0,0) (obviously V{0) = as well), then we get a Costas 
hypercube but we have no guarantee anymore that it has the Costas property; simulations show that 
sometimes it does. This is the equivalent of the Wo construction of Costas arrays by the addition of a 
"corner dot" to a W^i -constructed array [51[71IH]. 

• Welch arrays retain the Costas property when their columns get shifted circularly O [7]; this shift is 
expressed by the fixed parameter c in the definition of the permutation / shown earlier. As the extension 
of the Welch method formulated in Theorem |4] preserves this parameter c in the definition of /, we see 
that the hypercubes and the hyper-rectangles it produces have a similar periodicity: the fact that the 
family of dots at the locations {{i,f{i))\i £ [q — 1]} defines a Costas hyper-rectangle implies that the 
family of dots at the locations {(i, f{i ® k))\i G [g — 1]}, where i (B k :— 1 -I- [(i -I- fc — 1) mod {q — 1)], 
fc G [g — 1] — 1 fixed, also defines a Costas hyper-rectangle; similarly, the fact that the family of dots at 
the locations {{V (i) , f £ [q — 1]} defines a Costas hypercube implies that the family of dots at the 
locations {{V{i), f{i © k))\i G [g — 1]}, fc G [g — 1] — 1 fixed, also defines a Costas hypercube. 

• Just like we used V to obtain a Costas hypercube out of the originally constructed hyper-rectangle, we can 
use to obtain a permutation array of order g— 1, with dots at the positions {(i, {f(i)))\i G [g— 1]}. 
The question arises naturally whether this is actually a Costas array; alas, simulations show that it isn't. 
The reason is that does not generally preserve the Costas property, as we already noted in Remark^ 

• The most important aspect of this method is that it builds hypercubes not derived by Costas squares or 
arrays, at least in an obvious way. At the risk of sounding overly optimistic, if a method that converts 
Costcis hypercubes into Costas arrays were available, it could potentially lead to novel Costas arrays when 
applied on these hypercubes. 

5.3 A further generalization 

The isomorphism between the 2 representations of F(p'") as a vector space we suggested in Section 15.21 is 
probably the most obvious one, but by no means the only one possible. In terms of basis correspondence, note 
that V{p,m) is equipped with the natural basis of polynomials Pi{x) — x^, i G [m] — 1, F(p)"' is equipped with 
the natural basis of m-tuples — (5i,i, . . . , (5i,m-i), * G [m] — 1, and J-{Pi) ~ Si, i G [m] — 1. Alternatively, we 
could have used a different isomorphism J^b, such that the polynomial basis {Pi\i G [m] — 1} gets mapped to 
the rows of an invertible matrix B with elements in F(p), and consequently: 

TBiam^xx"^^^ + am-ix^^'^ + . . . + ao) = (bm_i, 6^-2, . . . , foo), "i, fei G F(p), i G [m] — 1 

where 

(flm-l, a.m-2, . . . , do) = (&m-l, bm-2, . • • , bo) B <^ (flm-l, am-2, • ■ • , ao)B ^ = (6m-l, 6m-2, . • . , bo) 

Therefore, J"s(-) = J^{-)B-^. 
Since B is invertible, 

{am-i, am-2, . . . ,ao) = {bm-i,bm-2, . • • ,6o) =0 

Theorem 5 (Welch hyper-rectangles and hypercubes under arbitrary bases). Let p be a prime, m G N*, <? a 
primitive root o/F(g) where q = p™, c G [g — 1] — 1, and B an invertible matnx with elements in ¥{p), so that 
its rows define a basis for the vector space F(p)'" over ¥{p). We choose Vp{p,m) as our representation o/F(g). 
Then: 

• The function fs '■ 1] F*(g), where fnii) = T{g^~^^'^ mod P{x))B~^ , i G [g— 1], P{x) an irreducible 
polynomial over F(p) of degree m, is a permutation over F* (g) . 

• The hyper-rectangle in m + 1 dimensions with side length g — 1 in the first dimension and p in the others, 
whose dots lie at the positions with coordinates {(i,/s(i))|i G [g — 1]}, has the Costas property. 

• The hypercube in 2m dimensions with side length p, whose dots lie at the positions with coordinates 
{(y(i), /s(i))|z G [g — 1]}, has the Costas property; V is the familiar mapping from Corollary [1^ with 
n — p m the present case. 

Proof. 

• /b is a permutation over F*(g) because g is taken to be a primitive root of F(g) and J^b is a bijection. 
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That the family of vectors {{i, /s(j))ii G [9 — 1]} has the Costas property follows from an almost verbatim 
repetition of the argument presented in Theorem^] consider the 4 integers ii, 12, 13 — ii + k, and 14 = i2 + k 
such that ij G [g- 1], j G [4], fe G N, ii < 12, and let (ii -12, fsiii) - /B(i2)) = {is-U, /b(«3) - /bCm)) <^ 
/fl(ji) — /s(22) = fB{iz) — fsiii)- This can be written as: 

■^i3 (3''"'"^" mod P(2;)) - Tb mod P(a;)) = 

= J^fl mod P(2;)) - Tb jj^^j p(^)^ ^ 

J^('<7'i"'+" mod P(2;))b~' - ^^(£1'=-'+" mod P(s)) B"' = 



jr(^^n+fe-i+c inodP(x)jB"'-J^(^<7'^+*-'+" mod P(s)j B"' <^ 
f mod P{x)) - T ig^'^""' mod P(x)) = 



and the rest of the proof proceeds directly as in the proof of Theorem |4l 

□ 

Remark 7. 

• It is clear that the ability to change the basis while preserving the Costas property increases tremendously 
the number of possible Welch-constructed Costas hypercubes. 

• All these hypercubes are "almost" permutation hypercubes: the addition of one dot at the zero position 
vector turns them into permutation hypercubes. 

5.4 Normal bases of fields and "rotational" constructions 

It is well known that in every finite field F(g), q = p™, p prime, m G N* there exists an element b such that the 
elements If , i = 0, . . . ,m — 1 are linearly independent over F(p), thus forming a basis of F(p)"' [T^. In the 
context of our generalized Welch construction, these bases (known as normal bases) lead to nice "rotational" 
hypercubes. The reason is the equivalence: 

m — 1 m — 1 m— 1 

i = i=0 i = l 

whereby x has coefficients (a;o, . . . , Xm-i) over the normal basis iff x^ has coefficients {xm-\,XQ, . . . , Xm-2), that 
is essentially the same expansion only cyclically shifted to the right. 

5.5 Examples 

Example 6. Let p = 3 and m = 3, so that q = 27, choose P{x) = a;^ + 2x + 1 which is irreducible over F(3), 
and choose c = 1, g = x. The Costas hyper-rectangle and the Costas hypercube constructed by Theorem |4] are 
shown in Table [G] along with the corresponding Welch permutation discussed in Remark |6l which fails to have 
the Costas property. In this particular example, adding a corner dot at (0, 0) to the hypercube preserves the 
Costas property, thus yielding a permutation Costas hypercube. 

Example 7. Let p = 5 and m = 2, so that q = 25, P{x) = x^ + x-\-2 which is irreducible over F(5), and c — 1, 
g — 2x. Further, choose the array B of Theorem[S]to be 



B = 



3 1 
2 



^ B-^ = 



3 1 
2 



The Costas hyper-rectangle and the Costas hypercube constructed by Theorem [5] are shown in Table [T] along 
with the corresponding Welch permutation discussed in Remark [H] which fails to have the Costas property. 
In this particular example, adding a corner dot at (0, 0) to the hypercube preserves the Costas property, thus 
yielding a permutation Costas hypercube. 
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Table 6: The process of constructing a Welch hypercube in F(27): the Welch hyper-rectangle corresponding to 
5 = .X, c = 0, P{x) = a;"^ + 2a; + 1 (left), the corresponding Welch hypercube (center), and the corresponding Welch 
permutation (right) 
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Table 7: The process of constructing a Welch hypercube in F(25): the Welch hyper-rectangle corresponding to 
g = 2x, c — 0, P{x) ~ + X + 2, B chosen as in Example [7] (left), the corresponding Welch hypercube (center), 
and the corresponding Welch permutation (right) 
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Table 8: The process of constructing a Welch hypercube in F(27): the Welch hyper-rectangle corresponding to 
g = 2x^, c = 1, P{x) = x'^ + 2x + 1 (left), the same hyper-rectangle after changing the basis using B chosen as in 
Example [8] (center), and the corresponding Welch hypercube (right), where we rearranged the order of the rows to 
demonstrate clearly the rotational structure 



Example 8. Let p = 3 and m = 3, so that q — 27, choose P{x) = + 2x + 1 which is irreducible over F(3), 
and choose c = 1, g — 2x^. It follows that = 2x^ + 2x + 2, g^ = 2x^ + x + 2, and one can see immediately 
that the 3 vectors {g,g'^,g^) are linearly independent. The matrix B corresponding to this choice of a (normal) 
basis is: 
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The Costas hyper-rectangle and the Costas hypercube constructed by Theorem [5] are shown in Table [8l In this 
particular example, adding a corner dot at (0, 0) to the hyper-rectangle (either before or after the change of 
basis) does not preserve the Costas property, but adding it to the hypercube does, thus producing a permutation 
Costas hypercube. 



6 Density of Costas hypercubes 

What is the maximal number of dots we can pack into a hypercube while retaining the Costas property? Up to 
this point we have avoided on purpose any discussion on this subject, precisely because it is a very complicated 
one. Note that in the case of Costas arrays this issue is not resolved intrinsically, but rather through the 
application of an extra condition: we want the array to be a permutation array, hence there can be exactly one 
dot per row and column. More generally, in a permutation Costas hyper-rectangle there are as many dots as 
the square root of the total number of elements in the hyper-rectangle. 
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This restriction has nothing to do with the Costas property itself: if we hft it (and thus allow any number 
of dots in a given row or column, even no dots at all, that is) , then we can easily find (through simulation) that 
we can create Costas squares of side n with more than n dots, even for values of n where Costas arrays are not 
known, like n = 32 or 33. Similarly, we can create hypercubes with more dots than the permutation restriction 
would allow, as we show below. 

6.1 Monte Carlo construction of Costas hypercubes 

The determination of the maximum number of dots in a hypercube under the Costas property alone seems to 
be a difficult combinatorial problem, and we have no solution for it at this time. Even in one dimension, this 
is the well known problem of the determination of the maximum number of dots in a Golomb ruler of a given 
length [16], which still remains unsolved. The best we can do is to try to construct Costas hypercubes through 
Monte Carlo simulation, using, for example, the algorithm below: 

Algorithm 1. 

1. Set the size n of the hypercube's side and its dimension m; 

2. Build a list of all possible n"" coordinate vectors in a random order; 

3. Place a dot in the hypercube where the first vector in the list indicates; 

4. For each coordinate vector of the list, from the second on to the last, test whether placing a dot in the 
hypercube where the vector indicates violates the Costas property: if it does, remove it; if it does not, 
leave it there. 

We hope that by running this Monte Carlo algorithm enough times we can obtain a Costas hypercube with 
a density close to the optimal one. Stochastic algorithms on Costas arrays, however, applied on more or less 
similar problems, are known not to work very well, precisely because Costas arrays are extremely rare [15) . 

Table [9] (left) shows the maximum number of dots in Costas hypercubes, for various values of the side length 
n and the dimension m, that we were able to construct by brute force, using the algorithm described above. 
Some of the hypercubes for which n = q — 2, q power of a prime, have been constructed by "sieving" the Golomb 
generalization, for comparison purposes: namely, instead of starting with the full list of coordinate vectors in 
the algorithm we described, we start with the list of vectors satisfying the first Golomb generalization in Section 
14.41 An important observation here is that using exclusively the latter vectors certainly restricts the potential 
of the construction method. 

6.2 Connection with 2 dimensions 

We can apply Algorithm [1] to construct Costas squares. Such constructions are given in Table |9] (right), where 
we see that for small n (approximately n < 50) the algorithm produces Costas squares more densely packed 
than Costas arrays, but this seems to be no longer the case for n > 60. In any case, having produced a Costas 
square for a given n, we can compare it to each Costas array of side length n and find the maximum number 
of dots that lie in common positions. A systematic high count would show that these arrays are pretty close to 
Costas arrays, so that, by perturbing the positions of a few dots and possibly removing or adding some others, 
we could easily obtain a Costas array. Unfortunately, in the tests we ran the count always turned out quite low. 

Therefore, in general Costas squares do not help us determine Costas arrays: although they typically have 
more dots than a Costas array requires, so that it might be expected that by carefully removing some the 
remaining ones would define a Costas array, in practice they almost invariably contain blank rows and/or 
columns (without dots, that is), whence it follows that, in addition to dot removal, it is still necessary to move 
dots around a bit into the blank columns and rows, and the way to do that is not obvious. Another way to 
state this is that, if we divide the dots of the array into equivalence classes consisting of either the rows or the 
columns of the array, choosing one representative only from each equivalence class almost invariably leads to 
incomplete permutations that cannot be completed in an obvious way while retaining the Costas property. 

7 Summary, conclusion, and future directions 

We defined the multidimensional generalization of Costas arrays in several possible ways, by investigating what 
the multidimensional generalization of the Costas property should be. We adopted the point of view that the 
Costas property of a multidimensional binary sequence depends exclusively on its autocorrelation, and that the 
permutation structure (which we also suitably defined in higher dimensions) is just an extra condition imposed, 
not directly related to the Costas property itself. 
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Table 9: The number of dots in some Costas hypercubes constructed by means of Algorithm [T] (left), as well as in 
some Costas squares (right). The ticked simulations (left) are produced by running Algorithm [T] not on the full list 
of position vectors, but only on those that satisfy the Golomb generalization in Section [44l 



Hence, Costas arrays became special cases of sequences with the Costas property in 2 dimensions, which we 
named Costas rectangles and squares; while the former generalized naturally to Costas hyper-rectangles and 
hypercubes, we generalized the latter in 2 distinct ways: 

• Strict Costas hypercubes: here, no 2 dots are allowed to have corresponding coordinates with the same 
value. We proposed 2 construction methods for such hypercubes, but we saw that the requirements of the 
definition severely limit the number of possible dots, and thus this case tends not to be very interesting. 

• Permutation Costas hypercubes: here, the hypercube was viewed as the representation of a permutation 
between vectors instead of integers. 

For the case of permutation Costas hyper-rectangles and hypercubes of even dimension, we proposed a 
construction method that reshapes an existing Costas array of suitable order into the desired hyper-rectangle 
or hypercube. In the case of odd dimension, we proposed a related heuristic that does not always work, but 
even when it doesn't it usually produces a hyper-rectangle that very nearly has the Costas property. We also 
generalized the constructions by starting with Costas squares instead of arrays, and gave specific examples. 

We subsequently investigated the application of the Welch construction method on finite fields with a non- 
prime number of elements, and found out that, although the method fails to produce Costas arrays, it produces 
Costas hyper-rectangles and hypercubes in a natural way, and actually in very large numbers, due to the possi- 
bility of changing the basis of the representation while retaining the Costas property. Once more, we supplied 
specific examples of the construction. 

Finally, we investigated experimentally, through Monte Carlo simulations, the difficult question of the restric- 
tions that the Costas property alone imposes on the number of dots in a Costas hypercube, and we generated 
directly Costas hypercubes and Costas squares of various side lengths in several different dimensions. We ob- 
served that, although in small side lengths the simulations produced hypercubes more densely packed with dots 
than those produced by the construction methods, as the side length increased this ceased to be the case. 

There are still many possible directions for future research in Costas hypercubes. For example, 

1. What is the maximum number of dots that can be packed into a Costas hypercube or hyper-rectangle, 
given the number of dimensions and the side lengths? 

2. Heuristic [T] seems to be producing Costas hypercubes pretty often, although in many occasions it fails to 
do so. Can we provide a rigorous and simple sufficient condition for the resulting hypercube to have the 
Costas property? 

3. Is there a different construction method for Costas hypercubes? In particular, can the Welch and Golomb 
methods be generalized in a direct way (that is, without the intermediate step of the use of the mapping 
V defines in Theorem [T]) in 3 or more dimensions? More generally, can a construction method be found 
directly based on finite fields? 

4. Are there any engineering applications of Costas hypercubes, perhaps of a similar nature to the applications 
of Costas arrays? 
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5. Is there a method to convert a Costas hypercube into a Costas array? Such a method could potentially 
lead to the construction of new Costas arrays, because some of the construction methods we have proposed, 
such as the extended Welch method in Section [S] produce (permutation) Costas hypercubes not linked to 
any Costas array, at least in an obvious way. 
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